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A quantum memory for light is a key element for the 
realization of future quantum information networks 
Requirements for a good quantum memory are (i) versatil- 
ity (allowing a wide range of inputs) and (ii) true quantum 
coherence (preserving quantum information). Here we 
demonstrate such a quantum memory for states possess- 
ing Einstein-Podolsky-Rosen (EPR) entanglement. These 
multi-photon states are two-mode squeezed by 6.0 dB with 
a variable orientation of squeezing and displaced by a few 
vacuum units. This range encompasses typical input al- 
phabets for a continuous variable quantum information 
protocol. The memory consists of two cells, one for each 
mode, filled with cesium atoms at room temperature with 
a memory time of about 1msec. The preservation of quan- 
tum coherence is rigorously proven by showing that the 
experimental memory fidelity 0.52 ± 0.02 significantly ex- 
ceeds the benchmark of 0.45 for the best possible classical 
memory for a range of displacements. 

A sufficient condition for a memory to be genuinely quan- 
tum can be formulated via the fidelity between memory in- 
put and output. If this fidelity exceeds the benchmark de- 
termined by the best classical device then it can store and 
preserve entanglement and hence is a true quantum memory. 
Classical benchmark fidelities are difficult to calculate, and 
until recently they were known only for coherent states 101 
and homogeneously distributed pure states jst]. Storage of 
non-classical correlations between single photons ||6|48[] and 
a dual-path superposition state of a photon |@| have been 
experimentally demonstrated with the Raman and electro- 
magnetically induced transparency (EIT) approaches. Along- 
side these demonstrations for single photons, quantum mem- 
ory for multi-photon entangled states is an important and chal- 
lenging ingredient of quantum information networks. 

A particular class of multi-photon entangled states, namely 
EPR-type entangled two-mode squeezed states, plays a fun- 
damental role in quantum information processing with con- 
tinuous variables (cv) ifloifTlll . A quantum memory for entan- 
gled cv states is valuable for iterative continuous variable en- 
tanglement distillation 11211. c ontinuous variable cluster state 
quantum computation il3Lll4ll . communication/cryptography 
protocols involving several rounds ifisll . and quantum illumi- 
nation \l^ 1/7]. Temporal delay of EPR entangled states has 
been demonstrated in lITsll . Very recently, classical bench- 
marks for storin g a s que ezed vacuum state lfl9ll and displaced 



squeezed states 11201 12111 have been derived. EIT-based mem- 



ory for squeezed vacuum has been recently reported 11221 12311 



albeit with the fidelity below the classical benchmark. 

Here we report the realization of a quantum memory for 
a set of displaced two-mode squeezed states with an uncon- 
ditionally high fidelity that exceeds the classical benchmark 

We store a displaced EPR state of two modes of light a+ 
and a_ with the frequencies uj± ~ ojq ± lol, where cjq is 
the carrier frequency of light. The entanglement condition is 
Var(X+ + X_) + Var(P+ - P_) < 2 jjH where canoni- 
cal quadrature operators obey [X±, P±] ~ i. For a vacuum 
state Var(Xvac) = Var(Pvac) = 1/2. The EPR entangle- 
ment of the a+ and a_ modes is equivalent to simultaneous 
squeezing of the cos(a;/,t) mode xlc ~ {X+ + X-)/V2; 
PLc = {P+ + P-)/ \/2 and the corresponding sin(cjLt) mode. 
The alphabet of quantum states is generated by displacing the 
two-mode squeezed vacuum state by varying values (xlcs) 
and {pLc,s) and by varying the orientation of the squeezed 
quadrature between XL andp^. The displaced squeezed states 
are produced (Fig. [T^) using an optical parametric amplifier 
(OPA) 1251] with the bandwidth of 8.3MHz and two electro- 
optical modulators (EOMs)(see the Methods summary for de- 
tails). 

The two photonic modes are stored in two ensembles of ce- 
sium atoms contained in paraffin coated glass cells (Fig. [T^) 
with the ground state coherence time around 30msec U]. ojq is 
blue detuned by A = 855 MHz from the F = 4 ^ F' 5 of 
D2 transition (Fig.[T];). Atoms are placed in a magnetic field 
which leads to the precession of the ground state spins with 
the Larmor frequency ujl = 2tt ■ 322 kHz. This ensures that 
the atoms efficiently couple to the entangled uj± = ujq ± ujl 
modes of light. The two ensembles 1(2) are optically pumped 
in F = A, rap = ±4 states, respectively, which leads to 
the opposite orientation of their macroscopic spin components 

Jxl — Jx2 — Jx- 

The atomic memory is conveniently described by two 
sets c, s of non-local, i.e., joint for the two separate mem- 
ory cells, canonical atomic operators ll26ll xac = [J^i^ — 

Jlf )/^,,PAc = {Jlf + Jlf )lV2J^XAs = -{Jlf - 

2Jx where the super- 



J^f)/^/2X,PA. = 
script rot denotes spin operators in a frame rotating at . It 
can be shown that the cosine/sine light mode couples only to 
the atomic c/s mode, respectively. As a consequence, in the 
protocol described below the upper(lower) entangled sideband 
mode of light is stored in the 1(2) memory cell, respectively. 
Since the equations describing the interaction are the same, 
we omit the indices c, s from now on. 
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Light sent from the sender station to the receiver (memory) 
station consists of quantum x-polarized modes and a strong y- 
polarized part which serves as the driving field for interaction 
with atoms and as the local oscillator (LO) for the subsequent 
homodyne measurement (Fig. [T]i. The interaction of light and 
a gas of spin polarized atoms under our experimental condi- 
tions can be described by the following input-output equations 

(1) 

where the coupling constant k is a function of light intensity, 
density of atoms and interaction time, and = 6.4 is a func- 
tion of the detuning only. In the limit n ^ Z these equations 
describe a swap of operators for light and atoms, i.e., a per- 
fect memory followed by squeezing by a factor Z^. However, 
in our experiment the swapping time is too long compared to 
the atomic decoherence time. We therefore follow another ap- 
proach which takes a much shorter time. 

The sequence of operations of the quantum memory proto- 
col is shown in Fig.[TJ). After the light-atoms interaction dur- 
ing the storage part of the protocol the output light operator 
x'^ is measured by the polarization homodyne detection and 
the result is fed back onto the p'^ with a gain g. This feedback 
into both c and s modes is achieved by applying two pulses of 
the magnetic field at the frequency oj^io the two cells. The 
resulting and p^ for the optimized g and k = 1 can be 
found from Eq. ([T]i 



a;i" = Y 1 - -^"^A + PL and p^" = -xl. (2) 

In the absence of decoherence the operator a;^ is perfectly 
mapped on the memory operator p^". The operator pL is 
stored in cc^" with the correct mean value (p^) = {x^) (since 
{xa) = 0) but with an additional noise due to xa- 

The ability to reproduce the correct mean values of the in- 
put state in the memory by adjusting g and k is a characteristic 
feature of our protocol which arguably makes it better suited 
for storage of multi-photon states compared to, e.g., EIT ap- 
proaches where such an adjustment is not possible. 

The additional noise of the initial state of atoms xa is sup- 
pressed by the factor — — 0.92 due to the swap in- 
teraction. To reduce it further we start the memory protocol 
with initializing the atomic memory state in a spin-squeezed 
state (SSS) with a squeezed xa- The SSS with variances 
Var (xa) = 0.43(3) and Var(p^) = 1.07(5) is generated 
as in [28] by the sequence shown as the preparation of initial 
state in Fig. [TJ) following optical pumping of atoms into the 
state with Var (xa) = 0.55(4). 

Before the input state of light undergoes various losses 
it is a 6dB squeezed state which we refer to as an "ini- 



tial pure state". In the photon number representation 
for the two modes the state is \^) = 0.8|0) + |G)_ + 
0.48|1) + |1)_ + 0.29|2)+|2)_ + 0.18|3) + |3)_ + .... We 
ran the memory protocol for 18 different states defined 
by = 0, 90 deg; [{xl) ; (pl>] - [0, 3.8, 7.6; 0, 3.8, 7.6]}, 
where (p is the phase of the squeezing, and (xl) and (p^) 
are the displacements of the initial state (see Table |l]i. In the 
absence of passive (reflection) losses for light and atomic 
decoherence we find from Eq. ^ the expected fidelity of 
0.95, 0.61 for the states squeezed with the = 0, = 90, 
respectively, with the mean fidelity of 0.78. 

In the experiment light is sent to the receiver's mem- 
ory station via a transmission channel with the transmis- 
sion coefficient r^ioss = .80(4) (which includes the OPA 
output coupling efficiency) resulting in the "memory input 
state" with Var (xl • cos{(f>) — pl ■ sin{(j))) = 0.20(2) and 
Va.r{xL ■ sin{(f)) +pl ■ cos{(f>)) ~ 1.68(9). The state of light 
is further attenuated by the factor T^cnt = .90(1) due to the en- 
trance (reflection) losses at the windows of the memory cells. 
Between the interaction and detection light experiences losses 
described by the detection efficiency ?7dct = .79(2) (see Meth- 
ods Summary). 

Following the storage time of 1msec we measure the atomic 
operators with the verifying probe pulse in a coherent state. 
The mean values and variances of the atomic operators x^ 
and p^" are summarized in Table|T](see the Methods summary 
for calibration of the atomic operators). From these values and 
the loss parameters, we can calculate the noise added during 
the storage process which is not accounted for by transmission 
and entrance losses. We find that the memory adds 0.47(6) 
to Var(x|") and 0.38(11) to Var(p|"), whereas for the ideal 
memory, according to Eq. (|2]), we expect the additional noise 
to be 0.36(5)(due to the finite squeezing of the initial atomic 
operator xa) and 0, for the two quadratures respectively. This 
added noise can be due to atomic decoherence, uncanceled 
noise from the initial anti-squeezed pA quadrature, and tech- 
nical noise from the EOMs. 

The overlap integrals between the stored states and the ini- 
tial pure states are given in Table U The average fidelities cal- 
culated from the overlap values for square input distributions 
with the size dmax = 0, 3.8 and 7.6 are plotted in Fig.|2] The 
choice of the interaction strength k = 1 minimizes the added 
noise but leads to the mismatch between the mean values of 
the stored atomic state and of the initial pure state of light by 
the factor -yy/JJ^^T^^^ioss = 0.85. This mismatch is the reason 
for the reduction of the experimental fidelity for states with 
larger displacements. 

The experimental fidelity is compared with the best classi- 
cal memory fidelity which is calculated f2B\ from the over- 
lap of the initial pure state with the state stored in the clas- 
sical memory positioned in place of the quantum memory. 
The input distribution used in the calculations is a square 
{K^l)! , \ {pl) \ < rfmax} with all input states with the mean 
values within the square having equal probability and other 
states having probability 0. The squeezing of the input states 
is fixed to the experimental value, and all phases of squeezing 



3 




Preparation of initial atomic state Storage Verification 



















1 


Optical 
Pumping 




Probe 1 




feed- 
back 




Probe 2 

/W 




feed- 
back 




optional 
Ti-pulse 




Probe 3 



no. 1 : Setup and pulse sequence, a. At the sender station two-mode entangled (squeezed) light is generated by the Optical parametric Am- 
plifier (OPA). A variable displacement of the state is achieved by injecting a coherent input into OPA modulated by electro-optical modulators 
(EOM). The output of the OPA is shaped by a chopper, and combined on a polarizing beamsplitter with the local oscillator (LO) beam, such 
that the squeezed light is only on during the second probe pulse. A beamshaper and a telescope create an expanded flattop intensity profile. The 
light is then send to the receiver (memory) consisting of two oppositely oriented ensembles of spin-polarized cesium vapour in paraffin coated 
cells and a homodyne detector. The detector signal is processed electronically and used as feedback onto the spins obtained via RF magnetic 
field pulses, b. Pulse sequence for the initiation of the memory, storage, and verification, c. Atomic level structure illustrating interaction of 
quantum modes with the memory. 



are allowed. The upper bound on the classical benchmark fi- 
delity is plotted in Fig. |2] The benchmark values have been 
obtained by first truncating the Hilbert space to a finite photon 
number and then solving the finite dimensional optimization 
employing semi-definite programming. The truncation of the 
Hilbert space is treated rigorously by upper bounding its effect 
via an error bound that is included in the fidelity. As a result, 
around the value dmax = 3.5 the theoretical calculated fidelity 
does not decrease further, but remains at a constant level due 
to the rapidly increasing truncation error (see the Supplemen- 
tary Material for the proof of the non-increasing character of 
the dependence of the benchmark on dmax)- 

For the input distributions around dmax = 3.8 (a vacuum 
unit of displacement is d = l/\/2), the measured fidelity is 
higher than the classical bound and thus the memory is a true 
quantum memory. 

Outperforming the classical benchmark means that our 
memory is capable of preserving EPR-entanglement in case 
when one of the two entangled modes is stored while the other 
is left propagating. Using experimentally obtained values of 
the added noise we evaluate the performance of our memory 



for the protocol where the upper sideband a+ mode is stored 
in one of the memory cells whereas the other EPR mode a_ 
is left as a propagating light mode. We find an EPR vari- 
ance between the stored mode and the propagating mode of 
1.52 (— 1.2dB below the separability criterion) which corre- 
sponds to the lower bound on the entanglement of formation 
of ^ l/7ebit ll29ll (see Supplementary Information for details 
of the calculation). Implementing this version of the memory 
would require spectral filtering of the a_ and d+ modes which 
can be accomplished by a narrow band optical cavity. 



In conclusion, we have experimentally demonstrated de- 
terministic quantum memory for multi-photon non-classical 
states. An obvious way to improve the fidelity is to reduce the 
reflection losses on the memory cell windows. Other improve- 
ments involve increasing the initial atomic spin squeezing and 
reduction of atomic decoherence which should allow for the 
storage protocol based on the swap interaction. 
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FIG. 2: Fidelities. The fidelity F calculated from experimental re- 
sults is shown as circles connected by lines. The theoretical bench- 
mark values are shown as squares. The horizontal axis is the size 
dmax of the input distribution. Note that one vacuum unit of dis- 
placement corresponds to dmax = l/%/2. 

Methods Summary 

Verification. By measuring the x'^ of the verification 
pulse the statistics of the atomic operator p^" is obtained. In 
another series we rotate into (and vice versa) using a 
magnetic 7r-pulse and obtain the statistics of x^ for the same 
input state of light. Since we assume Gaussian statistics, the 
mean values and the variances of x^ and are sufficient 
for a complete description of the atomic state. 

Calibrations. Before performing the storage, we calibrate 
the interaction strength k and the feedback gain g, such that 
the mean values of the light state inside the memory (i.e. af- 
ter the entrance loss) are transferred faithfully, k is calibrated 
by creating a mean value in the second probe pulse. 

The mean is stored in the atomic .t^, which is read out after 
the magnetic 7r-pulse with the third probe pulse. The mea- 
sured mean of the third pulse is then (a;'^"^) = (Pl"'^) 
from which K^is determined. Using similar methods we can 
calibrate the electronic feedback gain g. 

Generation of the displaced squeezed input states. The 
OPA which is pumped by the second harmonic of the master 
laser and generates the entangled squeezed vacuum states is 
seeded with a few /iW of the master laser light with the car- 
rier frequency ojo which is amplitude and phase modulated by 
two electro-optical modulators (EOMs) at a frequency of 322 
kHz, thus creating coherent states in the ± 322 kHz sidebands 
around (Fig. [T]i. With such a modulated seed, the output 
of the OPA is a displaced two-mode squeezed state. The out- 
put of the OPA is mixed at a polarizing beamsplitter with the 
strong local oscillator (driving) field from the master laser. 

Losses. In order to calculate the mean values and vari- 
ances of the stored state and the input states, we need to know 
the optical losses. We choose to divide the total losses j^tot 
into three parts, the channel propagation losses t^ioss from the 
OPA to the front of the memory cells (including the OPA out- 
put efficiency), the entrance loss r/ont and the detection losses 
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TABLE I: The three first columns display the mean displacements 
and the squeezing phase (0 = corresponds to being squeezed) 
of the initial pure light states. The next four columns display the 
mean values and variances of the atomic states after the storage. The 
last column displays the overlap between the initial pure light states 
and the stored atomic states. Vacuum state variances are 0.5. The 
uncertainties on the variances are calculated as the standard deviation 
of the variances within each subgroup of the input states. 



?7dct, such that 77tot=??ioss • '7cnt • ??dot (all the 77's are inten- 
sity transmission coefficients). From the measurement of the 
quadratures of the squeezed light (with variances 0.29(1) and 
1.34(6), we find the total losses jytot = 0.567(35). We mea- 
sure the transmission through the cells of 0.817(20), the trans- 
mission through the optics after the cells of 0.889(10) and es- 
timate the efficiency of the photodiodes to be 0.98(1). As- 
signing one half of the losses through the cells to the entrance 
losses and another half to the detection losses we find ?7ont = 
VCTST? = 0.90(1), ?7dct ^ 7(1817 • 0.889 • 0.98 = 0.79(2) 
andr^ioss = Vtot/ (?7cnt?ydct) = 0.80(4). 
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SUPPLEMENTARY INFORMATION 



Local/non-local operators and cos/sine vs sideband picture 



Memory added noise 

As discussed in the main text, in the experiment the final 
atomic states acquire extra noise from losses and decoher- 
ence not included in the basic theory. This extra noise is 
accounted for in two steps. In the first step the initial pure 
light state propagates through the transmission and memory 
entrance losses jyioss'yont (modelled by a beamsplitter) which 
add vacuum noise to the state of light. In the second step extra 
noise is added in the process of performing a unity-gain stor- 
age of the light state. We model this extra added noise with 
two operators Sx and Sp with zero mean values. The final 
atomic state is now (see also Eq. (2) in the main text) 



„fin 



Pa 



1 

Z2 



XA 



GpT" 



g^pT + Sx 



(3) 



where we defined the gain of the memory G = -/?/ioss?7cnt = 
0.85, which equals the ratio of the mean values of the stored 
atomic state and the intial pure state (this gain should not be 
mistaken with the electronic feedback gain), a;™"^ andpP'^ are 
vacuum operators with the zero mean and the variance 1/2. 
Taking into account that = 6.4 and the initial atomic noise 
is Var {xa) = 0.43 we can calculate the expected variances 
of the final atomic state from the above equation. 



,fin\ 



Var (x^' 



l_i_)Var(^^) + G2Var(pr°) 



Var [p^^) =G'Var (a;P" 



Var (5,), 
' 2 



+ Var(5p). (4) 



The variances of x^™'' 



and p^™'^ equal 1/ (2s)ors/2 depend- 
ing on the phase of the squeezing, where s = 4 in the exper- 
iment. By inserting these experimental parameters and the 
variances of the final atomic state, we can calculate the added 
noise. The results are given in Table Ullfor the two cases where 
we stored squeezed vacuum. 
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Below we will describe the quantum memory protocol in 
the language of local atomic operators and the upper and lower 
sideband modes instead of using the non-local atomic opera- 
tors and the cosine and sine light modes. The local atomic 
operators and the non-local operators are connected by 



1 



XAl [xAc +PAs 



XA2 



"7=2 ^""^ 



PAs 



PAl =-y| {PAc - XAs) , 

PA2 =-^{PAc + XAs) , (5) 



where the subscripts 1 and 2 refers to the first and second 
atomic ensemble, respectively. We can also find the relation 
between the upper and lower sideband modes of light and its 
cosine and sine mode 



XLc 



XLs 



1 

'72 

1 

'72 



{x+ +X-), 

[P- -P+) , 



PLC 
PLs 



{P++P-) , 
{x+ - X-) . 



(6) 



We write the input-output equations for the quantum memory 

as 



Xa 



-GplT + Ox, 



XAs 



^ pure 

--Gpl 



with the definitions 



PAc 
PAs 



GxlT 
GxlT 



Op 



p-l 



(7) 



1 



Ox=\\\- -^XA + - G^pT + Sx 



Op =v/l-G2 



Xl 



Sr, 



(8) 



Using Sx w 0.1 and Sp w 0.3, we find Var [Ox) ^ 0.44 and 
Var [Op) = 0.60. Equation (|7]i can be re-written by inserting 
Eqs. (|5]l and (|6]l, and we find 



Ox + Op 

Ox -Op 
V2 



Pai 



PA2 



GxP"' 



Op -Ox 
V2 ■ 
Op + Ox 



V2 



(9) 



We see that the upper sideband is stored in the first atomic 
ensemble and the lower sideband is stored in the second en- 
semble. 



TABLE IL The table shows the measured final atomic state variances 
Var [x^) and Var {pjC) and the added extra noise Var (Sx) and 
Var(Sp). 

The variances of Sx and Sp are quite small (less than one 
vacuum unit), and can be attributed to atomic decoherence and 
noise of the anti-squeezed light quadrature which can feed into 
the atomic p^" due to imperfect electronic feedback. 



Storage of one part of an EPR-entangled pair 

In the experiment, the sender prepares and sends a two- 
mode entangled state to the receiver, who then stores it in 
his quantum memory. Although the memory is a true quan- 
tum memory as proven by its ability to outperform the classi- 
cal benchmark, the noise added in the storage process to both 
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EPR modes leads to a separable state of the two memory cells. 
Note that for input states with the squeezing direction = 
displaced up to 3.8 (Table|I|, we are very close to having two 
displaced entangled atomic ensembles after the storage, since 
the parameter E describing EPR-entanglement 

E ^Var {xf, - x%) + Var {p^A + p%) 
=2.[Var(pi",)+Var(4",)] 
=2 -(0.52 + 0.52) = 2.08, (10) 

is only slightly above 2. 

Since outperforming the classical benchmark is sufficient 
to prove that the memory is capable of storing entanglement, 
we should be able to think of a modification of the experiment 
which can do exactly that. One example of such an experi- 
ment is the protocol, where Alice sends only one mode of the 
EPR-entangled pair to Bob for storage, and the other mode is 
sent to a third person, Charlie. In this case only one of the 
two entangled modes gets distorted by the memory. After the 
storage. Bob and Charlie perform a joint measurement to test 
whether there is entanglement between Bob's stored atomic 
state and Charlie's light state. In practice, Alice would have 
to separate the upper and lower sidebands, which could for 
instance be done using a cavity which transmits one sideband 
and reflects the other The initial entanglement of the upper 
and lower sidebands is characterized by 

E =Var (a;P"'" - a;P"") + Var (pP"" - pP"'=) = 2/s. 

(11) 

Alice then sends the EPR-entangled upper sideband (together 
with a lower sideband in the vacuum state) to Bob who stores 
the upper sideband in the first ensemble (and the vacuum in 
the lower sideband in the second ensemble). After the storage 
Bob and Charlie share the entanglement, since 

E =Var (xP"'" + p\\) + Var (pP"" - x\\) 
= ^ (1 + G) V ^ (1 - G) V Var (O,) + 
+Var(Op) = 1.5K 2, (12) 

which has been calculated using the experimentally obtained 
parameters of our memory from Table HH Eq. (|9]l and utilizing 
the fact that Var (a;P"" + x^_^'°) , Var (pP"'" + pP"'") = s. 
We conclude, based on the experimental performance of our 
quantum memory for the storage of both modes of the two 
mode entangled states, that if we instead stored only one part 
of the EPR-pair, one of the atomic memories would be entan- 
gled with the other part of the EPR-pair after the storage. 

Monotonicity of classical fidelity benchmark 

The following aims at proving that the classical benchmark 
is a non-increasing function of dmax- That is, the broader 




max 

FIG. 3: Fidelities. The fidelity F calculated from experimental re- 
sults is shown as circles connected by lines. The calculated bench- 
mark values are shown as squares and the monothonic benchmark 
values are shown as diamonds. The horizontal axis is the size dmax 
of the input distribution. Note that one vacuum unit of displacement 
corresponds to dmax = l/\/2. 

is the a priori distribution in phase space, the more difficult 
it gets for a classical scheme to achieve a certain fidelity. 
Here "classical scheme" refers either to entanglement break- 
ing quantum channels or to channels whose Choi matrix has a 
positive partial transpose. 

This supplement justifies the horizontal asymptote for the 
benchmark in Fig. 2 of the main text. Without this a rigorous 
treatment of the truncation errors (as in ref. Owari et al in the 
main text) would for large dmax lead to an increasing bound 
for the classical benchmark (see Fig.[3]l. 

Recall that our benchmark is based on the average fidelity 

F{T) := / d^qiO f ^ tr [T {Af^ipe,^)) pg.^] , (13) 

where ,t are the initial pure squeezed coherent states with 
mean ^ and orientation 9, J\f\ is an attenuation channel which 
decreases the intensity by a factor A, T is the channel corre- 
sponding to a hypothetical classical memory and g(^) is the 
probability density distribution of the input alphabet in phase 
space. The benchmark is then given by 

F := supF(T), 

T 

where the supremum runs over all T's corresponding to clas- 
sically possible schemes. 

We want to compare different values of F corresponding to 
different distributions q, so we regard F as a function of q and 
write Fq. Let us define 

T'ip) := £^ dfj p{v)W^T{W^^pWl^^)W,„ (14) 

where VF,, is the Weyl operator which displaces by rj in phase 
space and p{rj) is some probability density distribution. Note 
that T' characterizes a classical, albeit coarse-grained, scheme 
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if T does. In this case 



tain that 



i^, > fd^qiof f-tr[T'{Nx{pe,i))pe,i] (15) 

/do 
—di drj p{r])q{^)tT [T{Afx{p9,i+'n))p9.i+r,] 

= [d^q'iOf ^tr[T{N,{pe,^))pe,i] (16) 



(17) 



where we used that W^^Af\{p)W'^^^ = Wa (W,,pW,|) and 
we introduced the convolution 

Choosing T to be the optimal classical scheme for q' we ob- 



(18) 



under the assumption that q' is obtained from q by convolu- 
tion with another probability density p. If q and q', for in- 
stance, were two Gaussians, then Eq. (fTsT l holds whenever q' 
is broader than q since p can then be chosen to be a Gaussian 
whose variance is the difference between those of q' and q. 

Similar holds if, as in the experimental setup, q is a flat-top 
distributions on a square {~d, d] x {—d, d]. Taking the convo- 
lution with a discrete distribution = j ^iv ~ Vi)^ 
where the yy^'s are the four corners of the square [— d, d] x 
[— d, d], leads to a flat-top distribution q' on the square 
{—2d, 2d] X {—2d, 2d]. This justifies to use the benchmark 
computed for dmax = 3.8 again for dmax = 7.6 as done in 
Fig. 2 of the main text. 



